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Method of calculating hadron multi-point functions and disconnected quark loop contributions which are not 
readily accessible through conventional techniques is proposed. Results are reported for iv-n, tv-N and N-N 
scattering lengths and the flavor singlet-non singlet meson mass splitting estimated in quenched QCD. 



1. Introduction 

Numerical simulation of lattice QCD has been 
applied to an increasingly larger variety of strong 
interaction observables over the years. Yet these 
quantities all share the feature that their calcula- 
tion can be reduced to that of connected hadron 
2-point functions. A number of physically inter- 
esting quantities do not fall into this class: scat- 
tering amplitudes requiring hadron 4-point func- 
tions are a prime example. Amplitudes involving 
disconnected quark loops such as the flavor sin- 
glet meson propagator and tt-N sigma term repre- 
sent another important example. Technically the 
difficulty stems from the fact that calculation of 
these amplitudes requires quark propagators con- 
necting arbitrary pairs of space-time sites. With 
the conventional method of point source the nec- 
essary number of quark matrix inversions equals 
the space-time lattice volume, which would re- 
quire a prohibitively large amount of computer 
time. 

We found that the wall source technique Jl],^|], in 
particular without gauge fixing as were employed 
in the original proposals of extended sources |Q, 
could be effectively used to overcome the problem 
with a modest cost of computing power [[|. The 
method has been applied to calculate the full ir- 
7T, tt-N and N-N 4-point functions at vanishing 
relative momentum, from which we extracted the 



s-wave scattering lengths with the help of the re- 
lation between the two-particle energy in a finite 
box and the scattering lengthen . With a slight 
extension the method also allows an efficient cal- 
culation of the two quark loop contribution to the 
flavor singlet rf propagator. 

In this report we present the results together 
with some details of the technique. All of the sim- 
ulations have been made within quenched QCD 
at (3 = 5.7 mostly employing the lattice size of 
12 3 x 20. 

2. Calculational technique 

Consider the box diagram contributing to n- 
7r scattering with zero-momentum projected pion 
operators placed at the four time slices U(i = 
1, • ■ ■ , 4). To evaluate this amplitude on &n L 3 xT 
lattice, we calculate T quark propagators G t (n) 
with a wall source placed at the time slice t = 
1, ■ • • , T. With these propagators we form 



(^( G U^)G u (x 2 ,t 2 ) 
Gl(x 3 ,t s )G tl (x 3 ,t 3 ))), 



(1) 



*presented by Y. Kuramashi 



which equals the box amplitude except that the 
pion operators at the time slices t\ and have 
to be taken to be non-local in space without in- 
sertion of gauge link factors. 

Fixing gauge configurations to some gauge is 
usually employed to deal with the non-locality Q. 
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A problem with this method is that excited states 
such as p may contaminate signals for small time 
intervals. Alternatively one may take the gauge 
field average without gauge fixing since gauge vari- 
ant noise should cancel out in the average. One 
might worry that the noise overwhelms the signal 
in practice since for each time slice t\ and ti there 
are 0((L 3 ) 2 ) gauge dependent non-local terms 
relative to 0(L 3 ) local gauge invariant ones. It 
is our experience, however, that the noise level is 
controllably low at least for the pion. 

Calculating T quark propagators has another 
advantage. For two hadron operators placed at 
the same time slice, as is necessary if only the 
quark propagator for a wall source at t — is 
available, color Fierz rearrangement of quark lines 
takes place. This leads to a mixing among am- 
plitudes having different quark line topologies, 
which is numerically not straightforward to dis- 
entangle i|. Such a complication can be triv- 
ially avoided with the present method by placing 
hadron operators at different time slices. 

The present method can apply not only to a 
calculation of 4-point functions, but also to more 
general class of evaluation of n-point functions of 
hadrons. 

3. 7r-7r scattering lengths 

We have applied our method to calculate the 
1 = and 2 tt-tt scattering lengths for both 
Kogut-Susskind and Wilson quark actions, of 
which the Kogut-Susskind result has already been 
published jj). The procedure leading from tt-tt 
4-point functions to scattering lengths is as fol- 
lows. The energy shift of the two-pion state 
8E = E — 2™^ is extracted from the ratio of 4- 
to 2-point functions, 



R(t) 



< -rr(t + l)7r(f)7r(l)7r(0) > 



> 



< 7r(i)vr(0) 
Zexp(-SEt) 



< TT(t- 



(2) 



In practice we use a linear fit Z(l — SEt) in our 
quenched simulation since 0(t 2 ) terms are gener- 
ally not correct in the absence of dynamical quark 
loops. The fitted values can then be converted to 
s-wave scattering length ag using the Liischer's 




Figure 1. 1 = and 2 s-wave tt-tt scattering 
lengths Oq. Filled and open symbols denote Wil- 
son and Kogut-Susskind results. Triangles are for 
Coulomb gauge fixing. Squares (/3 = 5.7) and di- 
amonds (/3 = 6.0) for / = 2 are from Ref. [6]. 
Dotted lines indicate predictions of current alge- 
bra. 

relation (|]: 



E-2m w = 



^(l+ci— +c 2 (— ) )+0(L D )(3) 

nl-jx 1j J_j J_j 

with ci = -2.837297, c 2 = 6.375183. 

The two calculations respectively used 160 
(m q a = 0.Q1;KS) and 70 (K = Q.164; Wilson) 
12 3 x 20 configurations without gauge fixing. The 
1 = 2 amplitude has been studied in pioneering 
workj7j||. The calculation for / = is more diffi- 
cult since box and double annihilation diagrams 
with quarkless intermediate state, both requir- 
ing a full application of our method, contribute. 
We found that the double annihilation amplitude 
is very small, consistent with chiral perturbation 
theory and the OZI rule, and that the 1 = am- 
plitude exhibits a clear signal for attraction. 

Our results for the scattering lengths are sum- 
marized in Fig. |l| together with those of Ref. Q 
for 1 = 2. Triangles are the results obtained with 
Coulomb gauge fixing described below. We use 
TOjr and f n measured on the same set of configu- 
rations, corrected by the improved Z factor|| for 
/„■ for the Wilson case [Z = 1 for the KS case 
since the conserved current was employed). We 
observe an agreement of lattice results with cur- 
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rent algebra within 1-2 standard deviations for 
both 1 = and 2 channels. It is somewhat unex- 
pected that the agreement persists up to a quite 
heavy quark with m^/m p rj 0.7 — 0.8. 

The errors shown in Fig. |l| are statistical only. 
For the Kogut-Susskind action the lack of degen- 
eracy of pions in the Nambu-Goldstone and other 
channels invalidates the 0(L~ 5 ) term in (||)[[|. 
This leads to an uncertainty of 10% in the 1 = 
result ( that for the 1 = 2 result is less than 1%). 
Another source of systematic error is violation of 
scaling due to a fairly large lattice spacing of our 
simulation (a rj 0.2 fm at (3 = 5.7 ). The results 
of Ref. H obtained at j3 = 5.7 and 6.0 suggest, 
however, that this effect may be small, at least 
for 1 = 2. 

We have repeated the calculation with 
Coulomb gauge fixing employing 60 gauge con- 
figurations, and found the results for R(t) to be 
completely consistent with those obtained with- 
out gauge fixing on the same set of configura- 
tions. Comparing the two results, we found that 
the magnitude of errors is larger for the non-gauge 
fixed case. The amount of increase of errors, how- 
ever, is contained at the level of a factor of 1.5-2 
times those for the Coulomb gauge fixing, show- 
ing that gauge variant noise does not give rise to 
a serious problem for calculation of 7r-7r 4-point 
functions. The results for scattering length ob- 
tained with Coulomb gauge fixing are plotted in 
Fig.0. 

All the results in Fig. [l] are obtained in 
quenched simulations, and hence do not include 
effects of dynamical quarks. A potential prob- 
lem with the quenched calculation is that it could 
be affected by infrared singularities due to r\' 
loops suggested recently Since the appear- 
ance of such singularities is a general issue in 
quenched QCD, we shall discuss it in the subse- 
quent section, in connection with the U(l) prob- 
lem. The problem of this unphysical singularity 
anyway does not arise if calculations are made 
on full QCD gauge configurations [0, which is a 
straightforward task to be carried out. 

4. i]' meson mass in quenched QCD 

The problem of a large mass splitting between 
the flavor singlet rj and the pion octet is known 



as the U(l) problem. In the large N c expansion 
the mass splitting m\ = mz, — m 2 arises from 
iteration of virtual quark loops in the rj propa- 
gator, each loop giving a factor m 2 /(p 2 + m 2 ). 
In quenched QCD the series terminates at two 
quark loops. This implies that the magnitude of 
mo may be estimated through a comparison of 
the two quark loop amplitude having a double 
pole tuq/ (p 2 + m 2 ) 2 with the one quark loop am- 
plitude with a single pole l/(p 2 + JTj£); the ratio 
of the two amplitudes, each projected onto zero 
spatial momentum, is expected to behave as 



< rf(t)rf(fl) > 2 _ loop TO g 

= 7T~l ~ ~ 1. 



< tf(t)rf{0) > 1 _ loop 



2m T: 



(4) 



It has been suggested that the double pole in 
the rj propagator with a mass degenerate with 
that of ir leads to infrared divergences which are 
not suppressed by powers of the pion mass toward 
the chiral limit in quenched QCD||. To one- loop 
order in chiral perturbation theory, for instance, 
the pion mass receives a correction of the form, 



.1— loop\2 



1 



— — In- 

8N f n 2 f 2 A 



^ ■ (5) 



The problem for tt-tt scattering could be even 
more serious) 11 . The double annihilation dia- 
gram with two quark loops can be deformed into 
an rj' loop diagram with two insertions of the Trig 
vertex. Calculating the diagram in chiral per- 
turbation theory, one finds a divergent imaginary 
part at threshold in the s channel, and in the t 
channel a contribution of the form, 



SE = 



384tt 2 Nlm*f* L3 



(6) 



to the energy shift of the two pion system, which 
diverges quadratically as — » 0. Whether such 
terms affect results of quenched simulations de- 
pends on the magnitude of mo. A direct quenched 
estimate of mo through (^J) is also important in 
this respect. 

The two quark loop amplitude needed in (^) 
can be evaluated by the technique discussed 
above. We solve for the quark propagator with 
unit source at every space-time site without gauge 
fixing. (Note that this is crucial for the present 
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Figure 2. Ratio of two- and one-quark loop con- 
tribution to the rf propagator at K — 0.1665 and 
P = 5.7 on a 12 3 x 20 lattice. 

case.) With this quark propagator G(n,t) = 
X)(n" t „) G(n, t\ ri" , t") we form the expression, 



^Tr{G(n,0) 75 }Tr{Gt(n',i) 75 }. 



(7) 



This equals the two quark amplitude projected 
onto zero spatial momentum up to gauge-variant 
non-local terms which, however, cancel out in the 
ensemble average. A very nice feature of this 
technique is that it requires only a single quark 
matrix inversion for each gauge configuration. 

In Fig. ^| we plot the result for the ratio given 
in ([|) for Wilson quark action obtained with 160 
configurations at K — 0.1665. We observe a quite 
clear signal for a linear increase in t. We then 
extract the mo parameter by a fit of form (^) 
over the range 4 < t < 8. The results at two 
values of K together with a linear extrapolation 
in 1/ K to K c = 0.1694 are shown in Fig. || where 
conversion to physical units is made with a^ 1 = 
1.45(2) GeV determined from the p meson mass. 
We refer to earlier attempts jl2],[l3) that calculated 
the two quark loop amplitude. In particular the 
authors of Ref. |TJ| estimated mo = 530MeV at 
m^/mp = 0.71 and 330MeV at m^/rrip — 0.34 
with a -1 = 1.81GeV from a calculation with the 
rj' source and sink fixed at the spatial origin using 
10 configurations generated with a non-standard 
gauge action on an 8 3 x 16 lattice. 
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Figure 3. mo in physical units as a function of 
1/K - 1/K C . 

Our value m = 700(50)MeV at K c may be 
compared to the 'experimental' value of 850MeV 
deduced from the Witten-Veneziano formula [Q. 
Another comparison, more directly based on lat- 
tice QCD, is to invoke the original U(l) Ward 
identity relation m§ = 6%//^ with \ the topo- 
logical susceptibility. Using % = 4.28(33) x 10~ 4 
and /tt = 0.0616(40) at K c we obtained on the 
same set of configurations as for mo, we find 
m = 1190(90)MeV. 

Let us now consider the question of infrared 
singular terms induced by r( loops in quenched 
QCD. With our result for m and f n , we find 
that the coefficient of the logarithmic correction 
for the pion mass in ([5]) equals 0.076 at m^/nrip = 
0.61(if = 0.1665). Quenched pion mass data 
presently restricted to m^/m p ~ 0.5 do not show 
evidence of the logarithm with such a small co- 
efficient! 15 • The magnitude of (||) contributing 
to the 7r-7r scattering length is even smaller. For 
our Kogut-Susskind simulation, taking the Wil- 
son result mo a = 0.48 at K c as an indicative 
value and using measured and / w we find 
SE = 3.4 x 10~ 5 at m^/m.p — 0.32, which may be 
compared with the value 5E = 1.2(4.0) x 10~ 4 ex- 
tracted from the slope of the double annihilation 
diagram. We conclude that the possible failure 
of the quenched approximation due to rj' loops 
does not become manifest unless the simulation 
is made with a quark mass much smaller than the 
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value being taken in the current study. 

5. tt-N scattering lengths 

The wall source technique without gauge fixing 
does not yield good signals for the nucleon. We 
therefore used the Coulomb gauge fixing at the 
t = time slice for the nucleon source to calcu- 
late tt-N 4-point functions. In Fig. H we show 
our preliminary result for the ratio R (t) of tt-N 
amplitude divided by the propagators of tt and 
N obtained with 30 configurations with the Wil- 
son quark action at K = 0.164. We observe a 
good signal for I = 3/2. The data is much worse 
for / = 1/2, marginally indicative of a positive 
slope (i.e., attraction) expected from the experi- 
ment. Extracting the energy shift by a linear fit 
in t over the interval 4 < t < 8, we find for the 
scattering lengths 



7=3/2 




1=1/2 
J 



-0.95(13) 



T.6(0.7) 



-1 ] 



(8) 



(9) 



where jx^N = m^rriN / (m n + tun) and numbers 
in square brackets are predication of current alge- 
bra. Despite a heavy quark mass (m w /m p = 0.73 
at K — - 0.164) our results are consistent with cur- 
rent algebra both for I = 1/2 and 7 = 3/2 chan- 
nels. The error is large for I = 1/2, however. 

Clearly reducing the error is much desired for 
7=1/2. It is also interesting to repeat the 
calculation for smaller quark masses. However, 
both of these extensions are not easy in prac- 
tice. A simple argument suggests that the er- 
rors of R(t) grow exponentially SR(t) oc e at with 
a = tojv — 3/2??!^ or mjsr — l/2m„ depending 
on the topology of quark lines |l6[. Our numer- 
ical data indeed exhibits such dependence with 
the expected slope. This means that a reduction 
of errors at large values of t requires a substan- 
tial increase of statistics. The situation appears 
worse for lighter quarks since m^r remains non- 
zero while becomes smaller, leading to a larger 
value for the slope a. The only way to reduce er- 
rors would be to increase (3 with a simultaneous 
enlargement of the lattice size. In this regard the 
tt-N case turns out to be more difficult in nature 
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Figure 4. tt-N amplitudes for 7=1/2 and 3/2. 
Data are shifted by ±0.2 to avoid overlap. 



than the tt-tt case; for the latter the errors grow 
only modestly as 8R(t) oc e 2 ™"* even for the worst 
case of the double annihilation diagram || whose 
slope becomes smaller toward the chiral limit. 

6. N-N scattering lengths 

The N-N scattering lengths are experimentally 
known to be very large JL7J; 



ao( 3 Si) 



h20.1(4) fm, 
-5.432(5) fm. 



(10) 
(11) 



As opposed to the tt-tt and tt-N cases, the N- 
N amplitudes are not constrained by chiral sym- 
metry. The large values above are, therefore, a 
purely dynamical phenomenon whose successful 
derivation is a challenge posed to lattice QCD. 

We note that the negative sign in the 3 Si chan- 
nel is due to the presence of the deuteron bound 
state (Levinson's theorem). This brings in an ad- 
ditional complication in an attempt at a realistic 
calculation of the scattering lengths, since the en- 
ergy of the lowest scattering state orthogonal to 
the bound-state deuteron has to be computed to 
apply Luscher's relation to give the correct scat- 
tering lengths. 

Here we report on a less ambitious study 
employing quite a heavy Wilson quark with 
m^/irip = 0.86(-ftT = 0.16). For such a heavy 
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Figure 5. N-N amplitudes for 3 Si and 1 Sq chan- 
nels. Data are shifted by ±0.5 to avoid overlap. 



quark we may expect the deuteron to become 
unbound since the range of pion exchange is re- 
duced. In this case the scattering lengths for both 
3 5i and 1 5o channels can be extracted from the 
lowest N-N energy, and we expect both to be 
positive in sign and large. 

Anticipating a large scattering length we have 
made simulations on a 20 4 lattice with Wilson 
quark action using 30 gauge configurations with 
Coulomb gauge fixing. The result for the ratio 
R(t) is shown in Fig. [|. The positive slope corre- 
sponding to a negative energy shift for both chan- 
nels is consistent with our expectation above. Nu- 
merically, we find 



aoCSo) 
ao( 3 Si) 



1.0(3) fin, 
1.2(2) fin. 



(12) 
(13) 



for the scattering length in physical units us- 
ing a -1 = 1.45 GcV. Although these values are 
much smaller than the experimental values they 
are about a factor 3 — 4 larger than tt-tt and 
tt-N scattering lengths at similar quark masses, 
e.g., a I=a (tt-tt) = 0.421(23) fm and a I=1 / 2 {n- 
N) = 0.30(14) fin at K = 0.164. 

We consider that the present result is encour- 
aging; the most interesting problem left with us 
is to examine whether these scattering lengths in- 
crease rapidly as the quarks mass is reduced. 
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